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I.  INTRODUCTION 


A.  BACKGROUND 

A  pilot  must  be  furnished  with  flight  performance  information  for  his  aircraft 
that  will  permit  him  to  operate  the  aircraft  efficiently  at  different  airspeeds,  different 
altitudes,  different  values  of  gross  weight  and  different  external  loadings. 

This  information  is  collated  from  a  combination  of  theoretical  determinations 
and  flight  test  data  and  presented  to  the  pilot  in  a  flight  handbook  in  either  a  tabular 
format  or  in  a  series  of  graphs,  or  in  combinations  of  these  methods. 

B.  FUEL  FLOW  RATE 

Probably  the  most  important  parameter  in  optimal  performance  of  an  aircraft  is 
the  rate  at  which  fuel  is  consumed.  It  is  obvious  that  when  the  fuel  is  exhausted  the 
aircraft  will  no  longer  fly,  so  the  pilot  is  very  interested  in  the  rate  at  which  the  fuel  is 
being  consumed.  But  also  of  importance  is  the  fact  that  as  the  fuel  is  burned  the  weight 
of  the  aircraft  is  decreased  and  the  performance  of  the  aircraft  is  changed.  It  is 
therefore  quite  important  that  the  pilot  be  able  to  forecast  the  rate  of  fuel  consumption 
for  each  portion  of  each  flight. 

Fuel  flow  rate  for  a  turbo-jet  aircraft  is  a  direct  function  of  the  thrust  required 
for  the  specified  flight  condition,  and  in  equilibrium  flight  thrust  required  is  equal  to 
the  total  drag  of  the  aircraft. 

C.  TOTAL  DRAG 

The  total  drag  of  an  aircraft  in  equilibrium  flight  may  be  subdivided  into  induced 
drag  (the  drag  associated  with  the  generation  of  lift)  and  parasite  drag  (the  drag 
associated  with  the  resistance  of  the  fuselage  to  motion  through  the  air).  The  induced 
drag  is  a  function  of  the  altitude  of  the  aircraft  (air  density)  and  the  weight  of  the 
aircraft,  and  the  parasite  drag  is  also  a  function  of  the  altitude  as  well  as  being  a 
function  of  the  profile  drag  of  the  aircraft.  The  aircraft  drag,  which  varies  with  the 
configuration  and  external  loading  (weapons,  fuel  tanks,  et  cetera)  is  usually  expressed 
in  terms  of  a  Drag  Index  (DI),  with  the  zero  Drag  Index  usually  being  a  measure  of 
the  drag  of  the  aircraft  in  the  'clean'  configuration. 


D.  THRUST 


Although  it  is  the  thrust  of  the  entries  t.-.jt  p:  ...  ■ 
turbo-jet  aircraft,  as  long  as  there  is  sufficient  thrust  •  p-  . 

performance,  the  pilot  is  more  interested  lor  the  flight  plant.. r.g  .:.  v.c  : 

E.  FLIGHT  PLANNING 

To  plan  a  flight,  the  pilot  must  determine  what  will  be  the  loading  : 

(gross  weight)  and  what  effect  this  loading  will  have  on  the  drag  i  I)r  :g 
pilot  then  enters  the  flight  handbook  to  determine  the  fuel  How  rate  :  "  v 

and  altitude  at  which  he  desires  to  fly. 

If  the  gross  weight,  Drag  Index  and  airspeed  are  exactly  me  same  as 
in  the  handbook  tables  or  exactly  the  same  as  the  valued  on  the  graphical  laws.  tm 
determination  of  the  fuel  flow  rate  is  simply  a  matter  of  reading  (correctlx  i  the  value' 
from  the  information  present  in  the  handbook. 

If,  however,  any  of  the  input  parameters  are  not  identical  to  the  handbook  data 
points,  the  pilot  must  interpolate  between  the  values  presented.  If  weight.  Drag  Index, 
airspeed  and  altitude  are  all  off-presentation  points,  the  pilot  may  be  forced  to  do 
several  cross  interpolations.  This  may  even  involve  interpolations  between  values  listed 
on  different  pages.  Not  only  is  this  a  time  consuming  and  difficult  task,  but  the 
possibility  for  error  is  expanded  with  each  interpolation. 


II.  APPROACH  TO  THE  PROBLEM 


A.  STATEMENT  OF  THE  PROBLEM 

A  method  was  sought  that  would  permit  the  pilot  to  rapidly  obtain  the  fuel  flow 
rate  for  various  flight  combinations  of  gross  weight,  Drag  Index  and  airspeed.  And, 
regardless  of  the  ease  of  access  to  the  information  and  completeness  of  the  data 
obtained,  this  method  must  provide  fuel  flow  rate  at  the  same  degree  of  accuracy  as  the 
flight  handbook. 

In  addition,  it  was  planned  that  the  method  of  data  retrieval  should  be  as  simple 
as  possible  with  an  ultimate  goal  of  using  a  hand-held  computer  for  the  solutions. 

B.  BASE  LINE  APPROACH 

The  standard  method  of  using  the  flight  handbook  is  for  the  pilot  to  search  and 
follow  several  lines  and  tables  of  the  handbook  to  determine  the  required  fuel  flow  rate. 
This  requires  that  altitude,  gross  weight  and  Drag  Index  must  be  considered,  but  this 
method  is  too  long,  takes  too  much  effort  and  is  open  to  errors. 

The  procedure  can  be  simplified  by  computerizing  the  flight  manual.  Most 
attempts  at  this  method  have  involved  the  development  of  curve-fitting  routines  that 
will  supply  numerical  modeling  of  the  handbook  data.  Although  the  end  use  of  these 
equations  is  quite  rapid,  the  development  of  these  routines  is  a  long  and  arduous  task. 
Inasmuch  as  each  equation  usually  fits  only  a  small  part  of  the  data,  the  computer 
programs  for  such  a  method  are  usually  quite  lengthy. 

The  approach  undertaken  in  this  project  was  the  use  of  the  equations  for  an 
'IDEAL'  aircraft  to  develop  simplified  equations  that  can  be  used  in  a  short  computer 
program. 

Three  steps  were  undertaken  in  the  development  of  these  equations: 

a)  The  fuel  flow  rate  data  in  the  handbook  is  converted  to  thrust  required. 

b)  The  thrust  required  equations  are  adjusted  so  that  they  represent  the  actual 
aircraft. 

c)  The  corrected  thrust  required  equations  are  then  reconverted  to  flow  rate  as  the 
program  output. 
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C.  DETAILED  METHODS 
1.  Concept  of  Kj  and  K2 

Determinations  are  made  for  equilibrium,  unaccelerated  conditions  where  drag 
equals  thrust. 

The  basic  drag  equations  for  an  ideal  aircraft  are  used  [Ref.  l:p.  160] 

CD  =  CD0  +  CDj 

Cp  =  /Jt  e  AR 

D  =  (l/2)p  V2SCd 

All  of  the  principal  items  of  flight  performance  involve  steady  state  flight 
conditions  and  equilibrium  of  the  airplane.  For  the  airplane  to  remain  in  steady  level 
flight,  equilibrium  must  be  obtained  by  a  lift  equal  to  the  airplane  weight.  And  the 
airplane  drag  equals  the  thrust  required  to  maintain  steady  level  flight. 

The  total  drag  of  the  airplane  is  the  sum  of  the  parasite  and  induced  drag. 
Parasite  drag  is  the  sum  of  pressure  and  friction  drag  which  is  due  to  the  basic 
configuration  and,  as  defined,  is  independent  of  lift.  The  equation  of  parasite  drag  can 
be  expressed  as  a  function  of  squared  velocity. 

D0  «  (1  2)CDqPS  V2 

Let  Kj  =  (1/2)CDq  p  S 

D0=  K,V2 

Induced  drag  is  the  undesirable  but  unavoidable  consequence  of  the 
development  ol  lift. 

In  the  airfoil,  the  local  lift  vector  is  aligned  perpendicular  the  the  local  relative 
wind  and  hence  is  inclined  behind  the  vertical  by  the  angle  a  .  Consequently,  there  is  a 
component  of  the  local  lift  vector  in  the  direction  of  Vqq  ;  that  is,  there  is  a  drag 
created  by  the  presence  of  downwash. 

This  drag  is  defined  as  induced  drag,  denoted  by  Dj  [Ref.  l:pp.  152-155] 

Dj  =  (l/2)p  S  V2  CD 

Cp.  =  CL2  (ti  e  AR) 

Di  =  (CL>2  /(2n  c  AR)P  S  y2 
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=  ((2W/p  S  V1)1 1  (2k  e  AR)|  p  S  Vz 


=  (4W2  p  S  V2)  /  (p2  S2  V4  n  e  AR) 
=  [(2W2)  i  (p  S  7T  e  AR)](1/V2) 


If  the  whole  item  of  right  side  except  the  squared  velocity  is  changed  to  some 
constant,  this  equation  can  be  expressed  function  of  squared  velocity. 


Let  K2  =  (2W2)  /  (p  S  n  e  AR) 


=  K2(l/V2) 


While  the  parasite  drag  predominates  at  high  speed,  induced  drag 
predominates  at  low  speed.  Figure  2.1  illustrates  the  variation  with  speed  of  the 
induced,  parasite,  and  total  drag  for  a  specific  airplane  configuration  in  steady  level 
flight. 

The  curve  of  drag  or  thrust  required  versus  velocity  shows  the  variation  of 
induced,  parasite,  and  total  drag.  Induced  drag  predominates  at  low  speeds  when  the 
airplane  is  operated  at  maximum  lift-drag  ratio,  (L/D)max  the  total  drag  is  at  a 
minimum  and  the  induced  and  parasite  drags  are  equal. 

The  effect  of  a  change  in  weight  on  the  thrust  required  is  illustrated  by  Figure 
2.2.  The  changes  in  the  drag  curves  with  variations  of  airplane  weight,  configuration, 
and  altitude  furnish  insight  for  the  variation  of  range,  endurance,  climb  performance, 
etc.,  with  these  same  items. 

First,  the  primary  effect  of  a  weight  change  is  a  change  in  the  induced  drag  at 
any  given  speed.  Thus,  the  greatest  changes  in  the  curves  of  thrust  and  power  required 
will  take  place  in  the  range  of  low  speed  flight  where  the  induced  effects  predominate. 
The  changes  in  thrust  and  power  required  in  the  range  of  high  speed  flight  are 
relatively  slight  because  parasite  effects  predominate  at  high  speed.  The  induced  effects 
at  high  speed  are  relatively  small  and  changes  in  these  items  produce  a  small  effect  on 
the  total  thrust  required.  [Ref.  2:p.  99J 

Second,  Figure  2.3  illustrates  the  effect  on  the  curve  of  thrust  required  of  a 
change  in  the  equivalent  parasite  area  f  of  the  configuration.  Since  parasite  drag 
predominates  in  the  region  of  high  flight  speed,  a  change  in  f  will  produce  the  greatest 
change  in  thrust  and  power  required  at  high  speed.  Since  parasite  drag  is  relatively 
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Figure  2.2  EfTect  of  Weight  on  Thrust  and  Power  Required. 
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small  in  the  region  of  low  speed  flight,  a  change  in  f  will  produce  relatively  small 
changes  in  thrust  required  at  low  speeds.  The  principal  effect  of  a  change  in  equivalent 
parasite  area  of  the  configuration  is  to  change  the  parasite  drag  at  any  given  airspeed. 
[Ref.  2:p.  101] 

Third,  a  change  in  altitude  can  produce  significant  changes  in  the  curves  of 
thrust  required.  The  effect  of  altitude  on  these  curves  provide  a  great  part  c  '  the 
explanation  of  the  effect  of  altitude  on  range  and  endurance.  [Ref.  2:p.  101] 

Figure  2.4  illustrates  the  effect  of  a  change  in  altitude  in  the  curves  of  thurst 
required  for  a  specific  airplane  configuration  and  gross  weight.  As  long  as 
compressibility  effects  are  negligible,  the  principal  effect  of  increased  altitude  on  the 
curve  of  thrust  required  is  that  specific  aerodynamic  conditions  occur  at  higher  true 
airspeeds. 

As  being  investigated  above,  the  total  drag,  same  as  thrust  required,  is  the 
function  of  squared  velocity  with  different  constant. 

Total  drag  equation  is 

D  -  Kj  V2  +  K2  /  V2 
where  Kj  =  Cp  p  S 

K2  -  (2  W2)  /  (p  S  it  e  AR) 

The  first  step  of  this  project  is  to  determine  the  Kj  and  K2  values  to  obtain 
the  thrust  required.  It  can  be  seen  that  the  Kj  is  the  function  of  altitude  and  drag 
index,  K2  is  the  function  of  gross  weight  and  altitude.  So  if  K  j  and  K2  values  can  be 
computed  with  any  flight  condition,  the  thrust  required  will  be  solved  at  any  airspeed. 

2.  Determining  the  Kj 

To  determine  the  Kj,  use  the  minimum  drag  concept. 

At  the  minimum  drag  point,  the  airplane  can  be  flown  maximum  endurance 
and  the  induced  drag  and  parasite  drags  are  the  same. 

Th*  first  step  for  determine  the  K  j  is  to  investigate  the  thrust  equation.  The 
thrust  equation  can  be  established  with  fuel  flow  for  each  airspeed  and  some  factor  C. 

THRUST  -  FUEL  FLOW  x  C 

Since  the  fuel  flow  of  the  turbojet  powered  airplane  is  proportional  to  the 
airspeed,  the  thrust  required  is  proportional  to  the  fuel  flow. 
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Figure  2.3  Effect  of  Equivalent  Parasite  .Vea.f,  Thrust  and  Power  Required, 


The  factor  C  is  the  value  of  the  military  power  devided  by  fuel  flow'  at  military 

thrust. 

C  =  MILITARY  THRUST  /  FUEL  FLOW  (  at  military  thrust  ) 

The  military  thrust  and  fuel  flow  information  for  any  condition  can  be 
obtained  from  NATOPS  FLIGHT  MANUAL.  With  this  factor  C,  the  thrust  required 
for  each  speed  can  be  computed  by  fuel  flow, came  from  NATOPS  FLIGHT 
MANUAL,  multiplies  factor  C.  Table  1  and  Figures  2.5  and  2.6  show  the  relationship 
and  the  result  between  thrust  and  velocity. 

TABLE  1 

THRUST  REQUIRED  AND  FUEL  FLOW  VERSUS  VELOCITY 
VELOCITY(KTAS)  FUEL  FLOW(LBS/H)  THRUST  REQUIRED(LBS) 


360 

7,725 

6,558 

400 

9,014 

7,652 

440 

10,623 

9,018 

480 

12,462 

10,579 

520 

14,652 

12,438 

560 

17,377 

14,752 

600 

21  393 

18,160 

MIL 

25,680 

21,800 

Example)  military  thrust 

;  21800  lb 

fuel  flow  (at  mil)  ;  25680  lbs/hr 
Factor  C  -  21800/25680 
=  0.84891 


Thus 

Thrust  required  =  0.84891  x  Fuel  Flow 

Maximum  endurance  speed  at  specific  gross  weight,  altitude,  and  drag  index 
can  be  figured  out  with  the  graph  of  maximum  endurance  mach  number  in  NATOPS 
MANUAL  Figure  2.7.  This  true  mach  number  needs  to  be  transfered  to  KTAS  to 
substitute  into  the  drag  equation. 

Next  step  is  to  compute  the  minimum  drag. 

Minimum  drag  can  be  computed  from  the  fuel  flow  at  maximum  endurance 
multiplied  by  the  factor  C.  The  maximum  endurance  fuel  flow  can  be  obtained  from 
the  NATOPS  MANUAL  Figure  2.8. 

And  so  the  minimum  drag  is  the  twice  of  the  parasite  drag. 

^total  =  2^i  =  2Do  =  ^min 
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iguic  2.6  1  Iirust  Required  for  Velocity. 


Thus  the  drag  equation  can  be  expressed. 


W2  =  D0  -  (I'2)CDoPSv2 

where  p  =  density  at  given  altiude 
V  =  velocity  at  minimum  drag 
S  =  wing  area  of  specific  airplane 

Cn  can  be  obtained  from  the  above  equation,  because  all  values  except  Cn  are 
^o  uo 

known,  and  the  drag  equation  is  a  function  of  squared  velocity  with  the  known  Cp,  . 

The  total  value  of  right  side  except  is  the  desired  constant  Kj. 

3.  Determining  the  K-> 

The  procedure  to  compute  K2  is  just  similiar  to  that  for  computing  Kj.  It  is 
convenient  to  use  the  same  minimum  drag  concept.  At  the  minimum  drag  point,  the 
minimum  thrust  required  is  twice  the  induced  drag  amount. 

Use  the  same  maximum  endurance  velocity  and  minimum  thrust  required. 

The  induced  drag  equation  is 


Di  *  Tmin  '  2 


-  (2W 2)/(n  e  AR  p  S)(l/V2) 


where  w  =  Gross  weight  of  specific  airplane 
s  =  Wing  area  of  specific  airplane 
p  =  Air  density  at  given  altitude 
v  =  Maximum  endurance  velocity 


Since  the  minimum  thrust  required  and  all  items  of  the  right  side  are  known 
value  except  1/Jt  e  AR,  the  value  of  I  jit  e  AR  can  be  solved  easily.  After  obtained  the 
1  /tc  e  AR  value,  substitute  this  value  into  the  above  induced  drag  equation,  and 
calculated  the  numerical  value  of  the  right  side  except  the  l/V^.  This  numerical  value  is 
just  the  desired  constant  K2  value. 

But  if  the  raw  Kj  and  K2  values  are  used  to  compute  the  thrust  required, 
there  will  exist  a  deviation  of  more  than  10  percent  at  high  speed.  Therefore  one  needs 
to  modify  the  Kj  and  K2  to  fit  the  curve  to  the  FLIGHT  MANUAL  data. 

The  fitting  method  used  in  this  project  is  the  trial  and  error  method.  The 
Table  2  and  Figure  2.9  illustrates  the  difference  between  the  modified  and  raw  value  of 
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4.  Relationship  between  Kj,  K2  and  Variation  Factors 

Next  consideration  is  to  make  the  equations  of  Kj  and  K2  functions  of  drag 
index,  altitude  and  gross  weight. 

Since  the  K  j  is  the  function  of  drag  index  and  air  density,  K  j  will  be  changed 
by  different  drag  index  and  different  altitude.  Thus  it  is  needed  to  develop  a 
relationship  of  Kj  with  the  drag  index  and  altitude. 

To  accomplish  this,  it  is  convenient  to  select  some  reference  altitude  and  gross 
weight,  because  the  relationship  of  K  j  and  drag  index  is  independent  of  altitude  and 
gross  weight.  Each  Kj  value  for  every  drag  index  can  be  calculated  with  the  previous 
procedure.  After  calculation  of  all  of  the  Kj  values,  a  graph  of  Kj  versus  drag  index  is 
plotted  to  find  the  relationship. 

In  this  project,  the  altitude  and  gross  weight  are  selected  to  standard  sea  level 
and  40,000  lbs.  The  final  values  and  relationships  are  shown  in  Table  3  and  Figure 
2.10. 


TABLE  3 

DRAG  INDEX  VS  K1 


Drag  Index 

0 

20 

40 

60 

80 

100 

120 


K1 

0.011534 

0.013335 

0.014807 

0.016780 

0.018480 

0.020527 

0.022361 


Fortunately,  as  can  be  seen  in  the  result  of  the  graph,  the  relationship  between 
Kj  and  the  drag  index  is  linear.  So  the  equation  Kj  as  a  function  of  drag  index  is 


Kj  =  (0.0115  +  0.0001  x  Drag  Index) 

In  this  equation,  the  most  important  thing  is  the  slope  of  the  graph,  because 
the  starting  value  of  Kj  was  modified  to  fit  the  published  data. 

The  Kj  equation  with  modified  starting  value  of  Kj  is 

Kj  =  (0.014  +  0.0001  x  Drag  Index) 

This  equation  was  verified  by  calculating  the  drag  at  each  speed.  The  results 
are  shown  in  Tables  4,  5,  6. 
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TABLE  4 

THRUST  REQUIRED  DEVIATION  AT  DI  20 


ALTITUDE 


FT 


GROSS  WEIGHT 


40,000  LBS 


D.I 


20 


MILITARY 

THRUST 

21 

L , 800  LB 

THRUST 

FUEL  FLOW  x  0.8489 

FUEL  FLOW 

25,680  LB/H 

AIR  SPEED 

(KTS) 

FUEL  FLOW 

(LBS/H) 

THRUST  (LB) 

HAND  BOOK 

THRUST  (LB) 

CALCULATION 

DEVIATION 

m 

360 

7,725 

6,558 

6,456 

-1.5 

400 

9,014 

7,652 

7,741 

1.2 

440 

10,623 

9,018 

9,198 

2.0 

430 

12,462 

10,579 

10,819 

2.3 

520 

14,652 

12,438 

12,600 

1.3 

560 

17,377 

14,752 

14,536 

-1.5 

600 

21,393 

18,160 

16,625 

-8.5 

MAXIMUM  l 

FUEL  FLOW 

5,600  (LBS/H) 

ENDURANCE 

THRUST  (min) 

4,754  (LB) 

AIR  SPEED 

248  (KTAS) 

VARIATION  OF  K 

0.014+0.0001x20  =  0.016 

K2 

2.0  x  108 

From  Tables  4,  5,  6,  the  following  can  be  observed. 

First,  the  thrust  required,  i.e.,  calculated  drag,  is  almost  same  as  published 
data  at  low  drag  itdex  range,  but  the  thrust  required  is  much  less  than  published  data 
at  high  drag  index  range. 
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TABLE  5 

THRUST  REQUIRED  DEVIATION  AT  DI  60 


Second,  in  all  drag  index  ranges,  the  thrust  required  is  much  less  than 
published  data  at  high  air  speed  range  (approximately  above  Mach  number  0.8). 

To  compensate  for  the  first  problem,  that  is  low  thrust  required  in  high  drag 
index  range,  it  is  needed  to  increase  the  Kj  value,  but  the  thrust  required  value  is 
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TABLE  6 

THRUST  REQUIRED  DEVIATION  AT  DI  SO 


ALTITUDE  0  FT  GROSS  WEIGHT  40,000  LBS  D.I  80 

THRUST  21,800  LB 

MILITARY - -  — - THRUST  FUEL  FLOW  x0.9165 

FUEL  FLOW  23,785  LB/H 

AIR  SPEED  FUEL  FLOW  THRUST  (LB)  THRUST  (LB)  DEVIATION 


(KTS) 


(LBS/H)  HAND  BOOK  CALCULATION 


10,104 

12,094 

14,487 

17,419 

21,477 


9,261 

11,085 

13,278 

15,965 

19,658 


8,674 

10,479 

12,512 

14,763 

17,228 


-  6.3 

-  5.5 

-  5.8 

-  7.5 

-12.6 


MAXIMUM  ENDURANCE 


FUEL  FLOW 

6,650 

(LBS/H) 

THRUST  (min) 

6,003 

(LB) 

AIR  SPEED 

238 

( KTAS) 

VARIATION  OF  K 


0.014+0.001x80 
2.0  x  108 


=  0.022 


slightly  larger  than  the  published  data.  The  way  to  compensate  in  both  the  low  and 
high  drag  index  range  is  to  increase  the  slope  of  the  K  j  equation  and  to  decrease  the 
starting  value. 
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The  modified  equation  of  Kj  is 

Kj  =  (0.0135  +  0.000132  x  Drag  Index) 

And  the  change  of  Kj  equation  can  be  illustrated  in  Figure  2.11. 

The  second  problem,  that  the  thrust  required  is  much  less  than  the  FLIGHT 
MANUAL  data  in  the  speed  range  higher  than  Mach  number  0.8,  can  be  illustrated  as 
the  airplane  reaches  very  high  flight  speeds,  the  drag  rises  in  a  very  rapid  fashion  due 
to  compressibility.  Since  the  generalized  equation  for  parasite  drag  does  not  account 
for  compressibility  effects,  the  actual  dr  g  rise  is  typified  by  the  dashed  line  in  Figure 
2.12. 

This  phenomenon  can  be  investigated  from  the  viewpoint  of  drag  coefficient. 
The  lift  and  drag  coefficients  vary  with  the  Mach  number,  however,  the  drag  coefficient 
for  a  fixed  value  of  lift  coefficient  does  not  change  until  the  drag-divergence  Mach 
number  is  reached,  at  which  time  the  drag  coefficient  increases  sharply.  If  wind-tunnel 
tests  of  a  model  of  the  airplane  to  be  analyzed  have  been  conducted,  curves  of  drag 
coefficient  against  Mach  number  for  various  values  of  lift  coefficient  will  be  available 
as  shown  in  Figure  2.13. 

Drag  data  above  the  divergence  Mach  number  cannot  be  calculated  with 
precision,  and  to  achieve  any  accuracy  at  all,  wind-tunnel  tests  or  other  empirical  data 
must  be  relied  on.  It  is  noted  that  even  wind-tunnel  data  may  not  be  reliable  in  the 
transonic  range,  and  therefore  performance  guarantees  in  the  transonic  range,  based 
purely  on  wind-tunnel  data  and/or  calculations,  cannot  be  considered  completely 
dependable. 

In  most  instances,  the  lift-drag  relation  may  be  expressed  in  an  analytic  form 

such  as 

C_d  =  Cq  +  C /  7t  e  A R 

in  which  case  performance  characteristics  may  be  analytically  determined  if  similiar 
analytical  expression  for  power  available  can  be  written. 

The  equation 

CD  =  CD0  +  CL2  1  n  e  AR 

applies  to  most  aircraft  through  the  Mach  number  range  of  0  to  about  4  and  a 
range  from  0  to  0.6,  for  both  trimmed  and  untrummed  conditions,  however,  Crj  and  e 
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Figure  2.12  Typical  Airplane  Drag  Curve. 

become  functions  of  Mach  number  above  Mach  numbers  of  about  M  =  0.8.  For 
supersonic  aircraft  the  value  of  CD  increases  by  a  factor  of  2  to  3  through  the 
transonic  region  and  thereafter  generally  exhibits  a  slight  decrease.  [Ref.  3:pp.  275-278] 
The  span  efficiency  factor  e  shows  a  steady  decrease  from  M  =  0.8  to  about  M  =  2.0 
and  thereafter  remains  almost  constant  at  about  50  to  60%  of  its  low  speed  value. 
From  M  =  0.8  to  M=  1.2  the  parasite  drag  curve  quite  closely  resembles  a  sin^  curve, 
and  for  approximation  analysis  may  be  represented  by  the  following  equation 


KS 


The  variation  of  e  may  normally  be  fitted  by  the  polynomial 


e  =  eQ  -  e2  (  M  -  0.8  )2  +  e3  (  M  -  0.8  )3 


where  eQ  =  low  speed  value  of  e 

e2’  e3  =  constants  which  depend  on  the  airplane. 


With  the  above  theory  and  parasite  drag  equation,  it  is  possible  to  make  the 
drag  equation  between  the  Mach  number  M.8  to  M  1.2 

2 


D0  =  (l/2)PSCDo  W 


at  above  Mach  number  M.8,  where  Mach  Drag  exist, 


Dq  =  (l/2)p  S  [CD  +  (CD  -  CD  )sin3  {(M-0.8)  0.8}  1  V2 
u  ^oo  om  oo 


Kj  =  (1/2)  p  S  CD 


D0  =  [K,  +  K,/Cd  (Cd  -Cd  )  sin3  {(M-0.8)/0.8}  ]  V2 
u  1  1  ^oo  uom  oo 


-  K,  [  1  +  (CD  -  CD  )/CD  sin3  {(M-0.8)/0.8)  ]  V2 
1  lom  ^oo  oo 


Since  Kj  and  K2  are  function  of  air  density,  it  must  be  considered  altitude 
correction  factor  of  the  changing  altitude.  [Ref.  4) 

First  consideration  is  the  lapse  rate  and  vertical  temperature  structure  for  the 
standard  atmosphere. 

The  adopted  value  for  the  lapse  rate(-dT/dH  =  a)  in  the  troposphere  is  the 
following  constant 

a)  In  terms  of  standard  geopotential  meters 

a  =  0.0065°C  m'1 

b)  In  terms  of  standard  geopotential  feet 

a  =  0.0035661 6° F  ft’1 

c)  In  terms  of  the  c.g.s.  unit  of  H 


(a) 


a  =  6.628155  x  10'^  °C  cm'2  sec'2 


Accordingly  in  the  troposphere 


T  =  T0  -  all 


(b) 


where  a,H  and  TQ  are  in  consistent  unit 


Second  consideration  is  the  relationship  between  pressure  and  vertical 
displacement. 


The  following  relationships  are  to  be  understood  as  expressed  in  any  system  of 
consistent  units. 

The  air  is  assumed  to  obey  the  perfect  gas  law,  which,  in  the  c.g.s  system  of 
units,  may  be  written  as 

p  =  (1/RXP/T)  (c) 

where 

p  density  of  air,  gm  cm'^ 
p  pressure,  dynes  cm 
T  absolute  temperature,0  K 

R  gas  constant  for  1  gram  of  dry  air,  ergs  gm'1  {°K}'1 

And  the  air  is  assumed  also  to  be  in  gydrostatic  equilibrium  and  to  satisfy  the 
differential  equation 

dp  =  -  g  dz  (d) 

where,  in  c.g.s  units, 

P  pressure,  dynes  cm'^ 
p  density  (specific  mass),  gm  cm'^ 

Z  vertical  distance,  cm 
g  gravitational  acceleration,  cm  sec"^ 

The  vertical  displacement  is  herein  expressed  in  units  of  geopotential. 
Geopotential  is  defined  in  differential  form  by  the  equation 

G  dH  =  g  dz  (e) 

where 

Z  altitude  measured  positively  upwards  at  a  point 
g  positive  (absolute  numerical)  value  of  the  acceleration 
due  to  gravity  at  the  point 
H  geopotential  at  the  point 

G  dimensional  constant,  the  amount  of  which  determines 
the  magnitude  of  the  unit  of  H  in  terms  of  length  and 
time.(the  dimension  of  G  are  in  units  of  gZ  per  unit 
of  H.) 

Substituting  equation  (b)  in  equation  (c),  equation  (e)  in  equation  (d),  and 
combining  the  results  gives 


dP  P  =  G/R  (-dH/T0  -  aH) 


I 


=  G/aR  d(T0  -  aH)'(T0  -  aH) 

Intergrating  the  right-hand  member  between  limits  of  0  and  H  gives 

In(  P/P0  )  =  G/aR  ln(  TQ  -  aH  TQ  ) 

Let  n  =  G/aR 

Then 

P/P0  =  (T0  -  aH/T0  )n  =  (  T/T0  )n  (f) 

From  equation  (0,  using  the  consistent  numerical  values  of  G,a,and  R 
n  =  5.2561  (dimensionless) 

The  Figure  2.14  shows  the  relationship  between  temperature,  pressure,  density 
and  altitude. 

With  several  equation  investigated  in  above,  the  useful  relationship  between 
temperature,  pressure  change  and  vertical  displacement. 

For  this  project,  the  useful  relationships  are  written  as 

T/T0  -  1  -  6.875  x  10’6  H 
P/P0  =  (  1  -  6.875  x  1(T6  H  )5'2561 
T  =  518.688  x(  1  -  6.875  x  10’6  H  ) 
a2  =  1.4  x  1714.87  x  T 

Finally  the  relationship  between  density  of  air  and  vertical  displacement  is 

P/P0  =  (P/RT)/(P0/R0T0) 

=  (P/P0)(T0/T) 

And  it  can  be  seen  by  computing  several  cases  that  the  effect  of  changing 
altitude  becomes  larger  than  normal  density  changing  ratio. 

Let  K9  =(p/p0)1-2 

Then,  Kj  is  proportional  to  the  K9  and  K2  is  inverse  proportional  to  the  K9  . 

As  was  shown  in  previous  equations,  K2  is  proportional  to  the  squared  gross 
weight,  and  the  effect  of  changing  gross  weight  becomes  larger  than  normal  gross 
weight  changing  ratio,  same  as  altitude  change.  Thus,  if  there  was  a  specific  reference 
value  of  K2  for  specific  gross  weight,  the  relationship  between  arbitrary  gross  weight 
and  reference  gross  weight  can  be  made  easily  as 
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Tmc  ‘  N  Icinpci.iiuic.  ricNsuic,  Density  against  (icupotcntul  II. 


k2  =  Wre(4- 

With  given  three  specific  conditions, i.c.  altitude,  gross  weight  and  drag  index, 
the  thrust  required  can  be  computed  for  any  airspeed  condition  with  equations  shown 
above.  Kj  and  K->  values  can  be  computed  with  specific  altitude,  gross  weight  and  drag 
index,  and  the  thrust  required  for  any  air  speed  can  be  computed  with  Kj  and  K-, 

The  final  goal  of  this  project  is  to  compute  the  fuel  flow  rate  for  some  given 
condition,  i.e.  specific  altitude,  gross  weight,  drag  index  and  air  speed.  The  fuel  flow  is 
thrust  required  devided  by  factor  C.  The  factor  C  is  the  value  that  nulitarv  power 
devided  by  fuel  flow  of  that  time,  as  was  shown  in  previous  step. 

C  =  MILITARY  THRUST  FUEL  FLOW  (at  nulitarv  thrust) 

5.  Factor  C 

To  compute  factor  C,  it  is  necessary  to  find  the  relationship  rrulitarv  thrust 
change  for  different  altitude  and  fuel  flow  change.  This  is  done  in  two  steps. 

a.  Computing  the  military  thrust 

Altitude  is  one  factor  which  strongly  affects  the  performance  of  the  turbojet 
engine.  An  increase  in  altitude  produces  a  decrease  in  density  and  pressure  and,  ij 
below  the  tropopause,  a  decrease  in  temperature.  If  a  typical  nonafterbuming  turbojet 
engine  is  operated  at  a  constant  RPM  and  true  airspeed,  the  variation  of  thrust  and 
specific  fuel  consumption  with  altitude  can  be  approximated  from  Figure  2  14  I  he 
variation  of  density  in  the  standard  atmosphere  is  shown  by  the  values  of  densitv  ratio 
at  various  altitudes.  Typical  values  of  the  density  ratio  at  specific  altitudes  are  as 
follows 


Altitude,  ft: 


Densitv  ratio 


Sea  level  .  1.0000 

5,000  .  0.861' 

10,000  .  0.7385 

22,000  .  0.44 'b 

35,000  .  0.3044 

40,000  .  0  2462 

50,000  .  0.1532 

If  the  fixed  geometry  engine  is  operated  at  a  constant  V  i  I  ASi  in  subst  nic 
flight  and  constant  N  (RPM)  the  inlet  vclocitv,  inlet  ram.  and  compressor  pressure 
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ratio  are  essentially  constant  with  altitude.  An  increase  in  altitude  then  causes  the 
engine  air  mass  How  to  decrease  in  a  manner  very  nearly  identical  to  the  altitude 
density  ratio.  Of  course,  this  decrease  in  mass  flow  will  produce  a  significant  effect  on 
the  output  thrust  of  the  engine.  Actually,  the  variation  of  thrust  with  altitude  is  not 
quite  as  severe  as  the  density  variation  because  favorable  decreases  in  temperature 
o«.cur.  1  he  decrease  in  inlet  air  temperature  will  provide  a  relatively  greater  combustion 
gas  energy  and  allow  a  greater  jet  velocity.  The  increase  in  jet  velocity  somewhat 
offsets  the  decrease  in  mass  flow.  Of  course,  an  increase  in  altitude  provides  lower 
temperatures  below  the  tropopause.  Above  the  tropopause,  no  further  favorable 
decrease  in  temperature  takes  place  so  a  more  rapid  variation  of  thurst  will  take  place. 
1  he  approximate  variation  of  thrust  with  altitude  is  represented  by  Figure  2.15  and  can 
be  computed  with  the  equation  as  follows 

UvaJ-  hsI*  «"2) 
where  6  =  PaU  Pssl 
“  Palt  Pssl 

Some  typical  values  at  specific  altitudes  are  as  follows 

Altitude,  ft:  Ratio  or  (Thrust  at  altitude/thrust  at  SSL) 

Sea  level  .  1.0000 

5.000  .  0.8960 

10,000  .  0.8003 

15,000  .  0.7114 

20,000  .  0.6296 

25,000  .  0.5544 

30,000  .  0.4855 

Since  the  change  in  density  with  altitude  is  quite  rapid  at  low  altitude,  turbojet 
takeoff  performance  wall  be  greatly  affected  at  high  altitude.  Also  note  that  the  thrust 
at  ft.  is  approximately  39  percent  of  the  sea  level  value. 

The  thrust  added  by  the  afterburner  of  a  turbojet  engine  is  not  affected  so 
great bv  altitude  as  the  basic  engine  thrust.  The  use  of  afterburner  mav  provide  a 
'r.ru''  increase  of  percent  at  low  altitude  or  as  much  as  100  percent  at  high  altitude. 
Kef  2  pn  I  FM2i; 
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When  the  inlet  ram  and  compressor  pressure  ratio  is  fixed,  the  principal 
factor  affecting  the  specific  fuel  consumption  is  the  inlet  air  temperature. When  the  inlet 
air  temperature  is  lowered,  a  given  heat  addition  can  provide  relatively  greater  changes 
in  pressure  or  volume.  As  a  result,  a  given  thrust  output  requires  less  fuel  flow  and  the 
specific  fuel  consumption  is  reduced.  While  the  effect  of  altitude  on  specific  fuel 
consumption  does  not  compare  with  the  effect  on  thrust  output,  the  variation  is  large 
enough  to  strongly  influence  a  typical  variation  of  specific  fuel  consumption  with 
altitude.  Generally,  the  specific  fuel  consumption  decreases  steadily  with  altitude  until 
the  tropopause  is  reached  and  the  specific  fuel  consumption  at  this  point  is 
approximately  80  percent  of  the  sea  level  value. 

Above  the  tropopause  the  temperature  is  constant  and  altitude  slightly 
above  the  tropopause  causes  no  further  decrease  in  specific  fuel  consumption.  Actually, 
altitudes  much  above  the  tropopause  bring  about  a  general  deterioration  of  overall 
engine  efficiency  and  the  specific  fuel  consumption  begins  an  increase  with  altitude. 
The  extreme  altitudes  above  the  tropopause  produce  low  combustion  chamber 
pressures,  low  compressor  Reynolds  Numbers,  low  fuel  flow,  etc,  which  are  not 
conducive  to  high  engine  efficiency. 

The  available  military  thrust  for  this  project  with  various  altitude  is  shown 
on  Table  7  and  Figure  2.15. 

TABLE  7 

MILITARY  THRUST  FOR  DIFFERENT  ALTITUDES 


Altitude  (ft) 

Sea  level 
5.000 
1 6,000 
15,000 
20,000 
25,000 
30,000 


Military  thrust  (lbs) 

21,800 

19,533 

17,446 

15,508 

13,725 

12,086 

10,584 


The  appropriate  equation  for  the  curve  is  as  follows 

THRUST  =  21,800  -  0.4  x  ALTITUDE 

b.  Computing  the  fuel  flow  at  military  thrust 

There  are  three  effective  factors  required  to  change  the  fuel  flow.  One  is 
altitude,  another  is  drag  index,  and  the  other  is  gross  weight.  To  get  an  appropriate 
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fuel  flow  equation  as  a  function  of  altitude,  drag  index  and  gross  weight,  find  the 
relationship  between  fuel  flow  and  each  effective  factor,  then  combine  the  relationship 
totally. 

(1)  Relationship  between  fuel  flow  and  altitude. 

According  to  the  NATOPS  FLIGHT  MANUAL  data,  the  fuel  flow  is 
decreasing  gradually  with  increasing  the  altitude  with  other  conditions  constant.  To 
find  a  relationship  between  fuel  flow  and  altitude,  specific  drag  index  and  gross  weight 
are  fixed,  to  select  the  available  data.  Table  8  and  Figure  2.16  show  the  relationship 
when  drag  index  and  gross  weight  are  fixed  to  zero  and  40,000  lbs. 

After  curve-fitting,  the  fuel  flow  equation  stated  as  follows 

Fj  =  26,000  -  0.64  x  Altitude  (at  below  8,000  ft) 

Fj  =  20,840  -  0.37  x  (Altitude  -  8,000  )  (at  above  8,000  ft) 
where 

Fj  =  Fuel  flow  considered  with  different  altitude 

(2)  Relationship  between  fuel  flow  and  drag  index. 

According  to  the  NATOPS  FLIGHT  MANUAL  data,  the  fuel  Row  is 
decreased  as  the  drag  index  increased.  As  was  considered  in  step  (1),  a  specific  altitude 
and  gross  weight  will  be  fixed  to  find  the  relationship  between  fuel  flow  and  drag  index. 

Table  9  and  Figure  2.17  show  the  relationship  when  altitude  and  gross 
weight  are  fixed  to  zero  and  40,000  lbs. 

To  find  the  relationship  between  fuel  flow  and  drag  index  with  different 
altitude,  the  computed  fuel  flow  from  step  (1)  must  be  used.  With  above  two  variable 
coefficients,  the  fuel  flow  equation  is 

F'2  =  Fj  -  16.65  x  Drag  Index  (at  below  Drag  Index  20) 

F2  =  (Fj  +  333.3)  -  22.1  x  Drag  Index  (at  above  Drag  Index  20) 

where 

Fj  =  Fuel  flow  came  from  step  (1) 

F2  =  Fuel  flow  considered  with  different  altitude  and 
drag  index 

(3)  Relationship  between  fuel  flow  and  gross  weight. 

According  to  the  NATOPS  FLIGHT  MANUAL,  the  fuel  How 
decreased  a  small  amount  with  increasing  gross  weight.  To  find  the  relationship 


TABLE  8 

FEEL  FLOW  FOR  DIFFERENT  ALTITUDES 
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Altitude  (ft) 

Sea  level 
4,000 
8,000 
l  i.ooo 


Fuel  Flow  (lbs/hr) 

26,000 

23,520 

20,840 

19,000 


TABLE  9 

FUEL  FLOW  FOR  DIFFERENT  DRAG  INDEXS 


Drag  Index 


Fuel  Flow 

26,000 
25  680 
25  110 
24,425 
23,785 
23,225 
22,745 


TABLE  10 

FUEL  FLOW  FOR  DIFFERENT  GROSS  WEIGHTS 


Grossweight  (lbs) 

35,000 

40,000 

45,000 

50,000 

55,000 

60,000 


Fuel  Flow  (lbs/hr) 

26,000 

26,000 

25,995 

25,990 

25,985 

25,980 


between  fuel  flow  and  gross  weight,  a  specific  altitude  and  drag  index  must  be  selected 
as  in  the  previous  steps.  Table  10  and  Figure  2.18  show  the  relationship  when  altitude 
and  drag  index  are  fixed  to  standard  sea  level  and  zero. 

It  can  be  seen  from  Figure  2.18  that  the  relationship  between  fuel  flow 
and  gross  weight  is  linear.  The  final  equation  of  fuel  flow  with  different  three  variables, 
i.e.,  altitude,  drag  index  and  gross  weight,  can  be  written  as  follows 

F3  =  F2  -0.002x(Gross  weight-40,000)  (at  heavier  than  40,0001bs) 

F3  =  F2  (at  lighter  than  40,000  lbs) 


where 


F2  =  Fuel  Flow  came  from  step  (2) 
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F3  =  Fuel  Flow  considered  with  difTerent  altitude, 
drag  index  and  gross  weight 

Now  the  military  thrust  and  fuel  flow  for  any  altitude,  drag  index  and 
gross  weight  can  be  computed  with  the  above  step  (1)  and  (2).  The  desired  factor  C  can 
be  computed  with  the  military  thrust  and  fuel  flow. 

Factor  C  =  Military'  thrust/  Fuel  Flow 

To  compute  the  fuel  flow  for  any  various  specific  condition  is  the  final 
goal  of  this  project.  The  fuel  flow  can  be  computed  with  thrust  required,  i.e.  total  drag, 
devided  by  factor  C.  The  total  drag  and  factor  C  had  been  computed  in  previous  step 
already. 

FUEL  FLOW  =  THRUST  REQUIRED  /  FACTOR  C 

With  the  above  equation,  it  is  possible  to  find  the  required  fuel  flow  for 
different  specific  conditions  without  interpolating  the  flight  manual  data. 
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III.  RESULT 


A.  BASE  LINE  CONDITION 

Given  condition 

Gross  weight  40,000  lbs 

Drag  Index  20 

Altitude  Standard  Sea  Level 


1 .  Converting  velocity  versus  fuel  flow  to  velocity  versus  drag 

Military  thrust  (at  sea  level)  21,800  lb 
Fuel  flow  (at  military  thrust)  25,680  lb/hr 

Thrust  =  Fuel  Flow  x  Factor  C 
C  =  Military  thrust  /  Fuel  Flow 


at  military  thrust 

C  =  21,800/25,680  =  0.84891 

The  thrust  required  for  each  velocity  can  be  computed  by  fuel  flow  times 
factor  C(0. 84891).  The  Table  11  shows  the  result. 

TABLE  11 

THRUST  REQUIRED  FOR  EACH  SPEED 


VELOCITY(KTAS) 

360 

400 

440 

480 

520 

560 

600 

MIL 


FUEL  FLOW(LBS/H)  THRUST  REQUIRED(LBS) 


7,725 

9,014 

10,623 

12,462 

14,652 

17,377 

21,393 

25,680 


6,558 

7,652 

9,018 

10,579 

12,438 

14,752 

18,160 

21,800 


2.  Computing  the  minimum  thrust  and  maximum  endurance  speed. 

Maximum  endurance  speed  can  be  obtained  from  maximum  endurance  Mach 
number  graph  in  NAPTOPS  FLIGHT  MANUAL  (see  Figure  2.7).  The  maximum 
endurance  Mach  number  is  M  0.378.  And  maximum  endurance  fuel  flow  can  be 


obtained  from  the  maximum  endurance  fuel  flow  graph  in  NATOPS  MANUAL  (see 
Figure  2.8). 

The  minimum  thrust  can  be  obtained  from  the  maximum  endurance  fuel  flow 
by  multiplying  with  the  factor  C  computed  in  step  1. 

maximum  endurance  fuel  flow  ;  5,600  lb/ hr 
Tmin  =  5’600  x  °-84891  =  4,753.9  lb 

3.  Computing  the  K  j 

A'  maximum  endurance  speed,  the  thrust  required  is  equal  to  the  total  drag, 
that  is  twice  the  induced  drag  or  parasite  drag,  i.e.  induced  drag  and  parasite  drag  are 
the  same  at  this  point. 

^min  -  ^total  -  2^i  ~  2^o 
Do  =  CDo(l/2)pV2A 

where  D0  =  4,753.9  /  2  =  2,376.9  lb 
p  =  0.0023769  lb  sec2/ft4 
V  =  0.378  x  1,116.4  ft/sec 
A  =  530  ft2 

CDo  =  2Dq  /  (pV2A) 

=  (4, 753. 9)7(0.0023769  x  422.I82  x  530) 

=  0.02117183 

D0  =  (0.021 17183)(1/2)(0.0023769)(530)V2 
=  0.01333  V2 
Thus  Kj  =  0.0133333 

4.  Computing  the  K2 

Use  the  same  concept  as  for  Kj 

^min  =  ^total  =  2Di  =  2Do 
D0  =  2W2/(pAV2)(l/JteAR) 

1/rteAR  =  2W2  /  (PAV2DJ) 

=  (0.0023769X530)(422. 18)2(2376.9)/(2  x  400002) 

=  0.16677964 


Dj  =  (2W  pA)(  1  rteARXl  V*) 

=  (2  x  40,000"  0.0023769  x  530)(0. 166"964)(  1  V2) 


K2  =  4.2365  x  108 

Dt0tai  =  Thrust  required  =  KjV2  ~  Ki  V‘ 

These  Kj  and  K2  values  are  modified  a  little  to  fit  the  curve  of  hand  book 
data  by  the  trial  and  error  method  introduced  in  Chapter  2.  Thus  the  thrust  required 
equation  can  be  written  as 

D  =  0.01614V2  +  2.0xl08  V2 
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But  if  there  were  no  consideration  about  the  Mach  effects,  the  drag  at  high 
Mach  number,  i.e.  greater  than  M  0.8,  will  be  much  less  than  the  hand  book  data.  Cq 
becomes  a  function  of  Mach  number  above  Mach  numbers  of  about  M  0.8.  In  this 
region,  the  parasite  drag  curve  quite  closely  resembles  a  sin'*  curve,  and  for 
approximation  analysis  may  be  represented  by  the  equation 

CD  =  cd  +  <CD  -  CD  )  sin3  {(M  -  0.8), 0. 8} 

^o  oo  om  ^ oo > 

Thus  the  drag  equation  between  the  Mach  number  M  0.8  and  M  1.2 
D0  =  1/2  PSCDqV2 

Dq  =  1/2  pS  [CD  +  (CD  -  CD  )sin3{(M  -  0.8J/0.8}]  V2 
u  oo  uom  oo 

=  l/2pSCD 

D_  =  [K,  +  K,/Cd  (Cd  -  CD  )  sin3((M  -  0.8)/0.8}]  V2 
u  1  1  ^oo  uom  ^oo 

=  K,[1+(Cd  -Cd  )/Cd  sin3{(M  -  0.8)/0. 8} ]  V2 
om  ^oo  oo 

where  Cn  -  Cn  =  0.55 
uom  ^oo 

Cn  =  0.018 
00 

D0  =  Kj[l  +  0.55/0.018  sin3  {(M  -  0.8)/0.8}]V2 
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5.  Converting  to  fuel  flow  from  drag 

Fuel  flow  for  each  velocity  can  he  computed  from  thrust  required,  i.e.  total 
drag,  devided  by  factor  C. 

C  =  MILITARY  THRUST  FUEL  FLOW  tat  military  thrust) 

Military  thrust  and  fuel  flow  can  be  computed  with  the  equation  that  found  in 
Chapter  2. 

Military'  thrust  =  21.800  -  0.4  x  Altitude 
=  21,800  -  0.4  x  0 
=  21,800 

Fuel  flow; 

Fj  =  26,000  -  0.64  x  Altitude 
=  26,000  -  0.64  x  0 
=  26,000 

^2  “  Fj  *  16.65  x  Drag  Index 
=  26,000  -  16.65  x  20 
=  25,667 

F3  =  F2  -  0.002  x  (Gross  weight  -  40,000) 

=  26,667  -  0.002  x  (40,000  -  40.000) 

=  25,667 

Thus  C  =  21,800/25,667 
=  0.84934 

Fuel  flow  =  Thrust  required  /  0.84934 

Table  12  and  Figure  3.1  show  the  results  of  the  calculation  and  deviation  from 
hand  book  data. 

B.  VARIATION  NUMBER  1  -  GROSS  WEIGHT 
Given  condition 

Gross  weight  50,000  lbs 
Drag  Index  20 

Altitude  Standard  Sea  Level 


ALTITUDE  I  O  *’T  GROSS  WEIGHT  I  40.000  LBS  I  DRAG  INDEX 


320.0  360.0  400.0  440.0  480.0  520.0  560.0  600.0 

VELOCITY  (KTAS)  BASE  LINE  CONDITION 


The  gross  weight  is  changed  only  from  the  oase  line  condition.  It  is  necessary  to 
consider  the  change  of  Kj  and  Kv 

K,  =  1/2  pCDqS 

As  you  can  see  from  the  above  equation,  Kj  is  independent  of  gross  weight. 
Thus  there  was  no  change  in  Kj. 

But  K->  is  a  function  of  gross  weight  ( as  discussed  in  chapter  3). 

K2  =  2W2/(p  S  n  e  AR) 


It  can  be  seen  that  the  K2  is  proportional  to  the  squared  gross  weight,  but  the 
effects  of  varing  gross  weight  become  more  largely. 


K->  ■  (50.000.' 40,000)' 

^50, 0001b  240, 0001b 

=  2.0xl08  x(1.25)2  7 
=  3.653  x  108 


D  =  0.0164V2  +  3.653xIOS/V2 


Military  thrust  =  21,800  lb 
Fuel  flow; 

F  j  =  26,000  -  0.64  x  Altitude 
=  26,000 

F2  =  Fj  -  16.65  x  Drag  index 
=  26,000  -  16.65  x  20 
=  25,667 

Fj  =  F2  -  0.002  x  (Gross  weight  -  40,000) 

=  25,667  -  0.002  x  (50,000  -  40,000) 

-  25,647 

C  -  21,800/25,647 
-  0.85 

Fuel  flow  *  Thrust  required  /  0.85 

Table  13  and  Figure  3.2  show  the  results  of  calculation  and  deviation  from  hand 
book  data. 
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C.  VARIATION  NUMBER  2  -  DRAG  INDEX 

Given  condition 

Gross  weight  40,000  lbs 
Drag  Index  40 

Altitude  Standard  Sea  Level 

If  the  drag  index  was  changed  from  base  line  condition,  what  will  be  effect  on  the 
result?  It  is  necessary  to  consider  the  Kj  and  K2  equation. 

K,  =  l/2pSCDo 

K2  =  2W2  /  (pSrceAR) 

As  can  be  seen,  K2  is  a  function  of  gross  weight,  altitude,  wing  area  and  aspect 
ratio,  not  of  drag  index,  thus  K2  was  unchanged.  But  Kj  is  a  function  of  altitude,  wing 
area  and  drag  index.  Thus  the  appropriate  Kj  value  can  be  computed  with  equation 

Ki  =  K,  x (K.  /K,  ) 

*DI40  *nnn  '  lrw/wv  lmw 


‘DI20 


DI40  *DI20 


Now  it  is  necessary  to  find  the  relationship  of  Kj  between  drag  index  20  and 
drag  index  40.  Furthermore,  if  the  relationship  of  Kj  between  any  drag  index  and 
reference  drag  index,  the  various  Kj  for  different  drag  index  can  be  computed.  To 
figure  out  the  relationship  of  K  j  for  each  drag  index,  let  us  compute  each  K  j  value. 
The  computing  method  is  the  same  as  introduced  in  the  base  line  condition  step. 

Tmin  =  2Do 
D0  =  1/2  PSCDqV2 

CDq  =  2D0/(pSV2) 

Kj  =  I/2pSCDq 

For  example 
when  DI  ■  0 

D0  =  T^/2  =  4,443.8/2  =2,221.9 

CD  =  2(2,221.9)/(0.0023769x530x438.92) 

°  =0.018132 


K,  =  l/2pSCD 


^ <■  :  v.  v.vr;  v;.- ;  ,»v  .■ 


v  .'-j.  *  SVv'  •'w 


=  1/2  (0.0023769  x  530  x  0.018132) 
=  0.011534 


V 
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I 
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Table  14  and  Figure  3.3  show  the  result  and  the  relationship  of  Kj  for  different 
drag  index. 

TABLE  14 

K1  VALUE  VERSUS  DRAG  INDEX 


DI 

c 

MAXIMUM  ENDURANCE 

THRUST (LB) 

CDo 

FUEL  FLOW 

KTAS 

MINIMUM 

K1 

0 

0.8385 

5,300 

260 

4,444 

0.018 

0.012 

20 

0.8489 

5,600 

250 

4,754 

0.021 

0.013 

40 

0.8682 

5,920 

247 

5,140 

0.024 

0.015 

60 

0.8925 

6,250 

241 

5,578 

0.027 

0.017 

80 

0.9165 

6,520 

238 

5,976 

0.029 

0.018 

100 

0.9386 

6,800 

233 

6,383 

0.033 

0.021 

120 

0.9584 

7,050 

230 

6,757 

« 

0.036 

0.022 

The  K  |  of  the  base  line  condition  was  modified  to  fit  the  drag  index.  As  you  can 
see  in  Figure  3.3  and  the  investigation  in  Chapter  2,  the  relationship  between  Kj  and 
drag  index  is  linear  and  the  equation  of  K  j ,  function  of  drag  index  is 

Kj  -  0.0135  +  0.000132  x  Drag  Index 
Thus 

K,  =  0.0135  +  0.000132  x40 
DI40 

=  0.01878 

Thrust  required  =  0.01878V2  +  2xl08/V2 
Military  thrust  =  21,800  lb 


mm 


Fuel  flow; 


Fj  =  26,000 


K- 


* 


F2  =  (Fj  +  333.3)  -  22.1  x  Drag  Index 
=  (26,000  +  333.3)-  22.1  x  40 
=  25,449.3 


C  = 


F3  =  F2  -  0.002  x  (40,000  -  40,000) 
=  25,449.3 

21,800/25,449.3 

0.8566 


Fuel  flow  =  Thrust  required/0.8566 


Table  15  and  Figure  3.4  show  the  results  of  calculation  and  deviation  from  hand 
book  data. 

D.  VARIATION  NUMBER  3  -  ALTITUDE 

Given  condition 

Gross  weight  40,000  lbs 
Drag  Index  20 

Altitude  4,000  ft 

If  the  altitude  will  be  changed  only  from  the  base  line  condition,  how  will  it  affect 
the  Kj  and  K2? 

The  equations  of  Kj  and  K2  are 
Kj  =  1/2pSCDq 
K2  -  2W2/(pSrceAR) 

As  can  be  seen  in  the  above  equation,  the  Kj  and  K2  both  are  functions  of  the 
air  density,  and  air  density  is  a  function  of  altitude.  Thus  the  Kj  and  K2  must  be 
changed  with  varying  altitude. 

The  useful  relationships  came  from  the  ICAO  Report  written  as 


p  =  P/RT 

T/Tq  =  1  -  6.875x10'6H 
P/P0  =  (1  -  6.875x10’6H)5-2561 


hi 


65 


MAXIMUM  ENDURANCE  THRUST  (MINIMUM)  5,920  x  0.86818  =  5,140  (LB) 


INDEX  VARIATION 


520.0  500.0 

Ac;  INDEX  VARIATION 


T  =  518.688  x  ( 1  -  6.875xlO‘6H) 

And  the  relationship  between  density  of  air  and  vertical  displacement  is 

p/p0  =  (P/RTXP0/RT0) 

=  (P  P0XT0/T) 

Since  the  effects  of  changing  altitude  become  larger,  it  is  necessary  to  modify  the 
equation. 

Let  K9  =  (p/p0)1'2 

Then  Kj  is  proportional  to  the  K9  and  K->  is  inverse  proportional  to  the 
Thus  in  this  case 

PPo  =  (PP0KT0T) 

=  (  1  -  6.875xl0'6x4.000)5  2561  (1  -  6.8”’5xlO'6  x  4,000) 

=  0.88809 

K  ■  =  K,  x  (0.88809)1'2 
ssl 

=  .0139976 

K-,  =  K->  (0.88809) 1  2 

2  ■‘■ssl 

=  2.306xl08 

Thrust  required  =  0.0|  399"r6V'  *  2  30bxl<>8  V*- 


Ihe  mihtarv  thrust  will  be  reduced  as  discussed  in  <  hapter  2 


Tava,l  =  1  ssi  6  ^ 


,«1  2i 


Pj.oooft  Pssl 

-  1 .82"  W  2.1  I*  22 

-  o  H6 tftS'H 


®  "  P4  1  P vs ; 

m  I  I  1  X  l  2  1  1  *  1  11  *  1  2  '  ) 

-  ■  1  V  X  Si  >  .  ‘ 


.1  a 


• »  t 


«  •  V  '  V  4 


f  uel  How; 


Fj  =  26,000  -  0.64  x  Altitude 
*  26.000  -  0.64  x  4,000 
=  23.440 

F2  =  ^  |  *  16.65  x  Drag  Index 
=  23.440  -  16.65  x  20 
=  23,107 

F3  =  F2  -  0.002  x  (Gross  weight  -  40,000) 

=  23,107  -  0.002  x  (40,000  -  40,000) 

=  0.86463 

C  =  19.979  23,107 
-  0.86463 

Fuel  (low  =  Thrust  required'O. 86463 

Table  16  and  Figure  3.5  show  the  results  of  calculation  and  deviation  from  hand 
book  data. 

E.  V  ARIATION  NUMBER  4  -  WEIGHT,  DRAG  INDEX  AND  ALTITUDE 

Given  condition 

Gross  weight  50,000  lbs 
Drag  Index  40 

Altitude  8,000  ft 

If  :*M  conditions  are  changed  from  base  line  conditions,  it  is  neceassary  to 
consider  for  each  changed  condition  step  by  step,  same  as  previous  sections. 

At  first ,  think  about  gross  weight  change  rhe  changing  gross  weight  affects  the 
k->  oni\.  because  K|  is  independent  of  gross  weight. 

K2  -  2W2  1  pSneAR) 

\t  you  can  see  in  section  *  2 

k -*  »  K,  i  Viooo  40,000  r 

*  S< (  t XU  ilbs  *40.0O0|bs 

-  2  ox|os  i  1  2S>‘ 

-  ’ 


muon. 


Next,  think  about  the  drag  index  change.  If  the  drag  index  is  changed,  it  is 
necessary'  to  consider  K|  only,  because  the  is  independent  of  the  drag  index  change 
K  |  equation,  function  of  drag  index  is 

Kj  *  0.0135  +  0.000132  x  Drag  Index 

Thus 

K  |  -  0.0135  +  0.000132  x  40 

DI40 

-  0.01878 

And  finally,  think  about  the  altitude  change. 

As  you  can  see  in  section  3.4,  the  altitude  change  afreets  both  Kj  and  K2. 
because  Kj  and  K2  are  function  of  air  density.  The  equation  in  section  3.4  can  be  used 
in  this  step. 

P;P0  -  (P'RTXP0/RT0) 

"(P/PqXTq/T) 

Let  K9  -  (P  Pq)1-2 
Kj  -  f(P) 

K2  -  f  (1/p) 

Thus 

K1  *  KlD140x(P8,000ft  Pssl^  2 

-0.01878  x  (0.0018683,0.0023769)*  2 

=  0.01406755 

K2  ’  ^O.OOOIbs  '  (P8-0°Oft  P”'}' 2 

-  3.653xl08  (0  0018683  0.0023769)*  2 

•  2.9 157x1 0® 

Throat  required  -  001407V2  ♦  2  9|S-XI08  V2 

I  he  militarv  thrust  will  be  computed  a  -  discussed  in  thaler  2 
t  x  I  2 

1  x.nooft  “  1  ssl  ,d  a  1 

*  ”  Ps  nnoft  *\sl 


8,000ft 


=  1571.88  21 16.22 
=  0.7427',7 

0  *  P8, 000ft  Pssl 

-  0.001 8683  0. 0023769 
=  0.786024 

-*  21,800  x  (0742777  0.786024 1  2) 


=  18.2641b 


Fuel  flow; 

F  i  *  26,000  -  0.64  x  Altitude 

-  26.000  -  0.64x8,000 

-  20,880 

^2  "  (F|  +  333.3)  -  22.1  x  Drag  Index 

-  21.213.3  -  22.1  x  40 

*  20,329.3 

Fj  *  Fj  -  0.002  x  (Gross  weight  -  40,000) 

-  20.329.3  -  0.002  x  (50,000  -  40,000) 

-  20,309.3 

C  »  Military  thrust  '  Fuel  flow 
*  18,264  20,309.3 
-  0.8993 

Fuel  flow  -  Thrust  required/0.8993 

fable  17  and  Figure  3.6  show  the  results  of  calculation  and  deviation  from  hand 
book  data. 

F.  VARIATION  NUMBER  5  -  WEIGHT,  DRAG  INDEX  AND  ALTITUDE 

We  have  wen  how  the  fuel  flow  was  changed  with  various  conditions  Not  much 
time  is  needed  to  find  the  fuel  flow  for  conditions  which  can  be  read  directly  from  hand 
book,  but  most  of  the  conditions  do  not  correspond  to  the  published  conditions  which 
can  be  read  directh  from  the  hand  book  data  sheets  Thus  to  find  the  fuel  flow  tor 
specific  conditions  that  are  not  printed  in  the  hand  book,  one  must  spend  lots  of  time 
’  •  interpolate  each  item  I  he  specific  fuel  flow  will  now  be  computed  b\  both  methods 


*  9 


VAR  I A 


320.0  360.0  400.0  440.0  480.0  520.0  560.0 

VELOCITY  (KTAS)  ALL  CONDITION  VARIATION 


- 1 - 1 - 1 - : - - - 1 - ; - ; 

0#SI  0‘tt  OCT  OST  011  0*01  0‘6  0‘9  07. 

tOI.  (HH/Sai)  MOU  13HJ 


Given  condition 


Gross  weight  48,000  lbs 

Drag  Index  27 

Altitude  11,000  ft 

Velocity  455  KTS 

The  procedure  to  compute  the  Kj  and  K2  is  the  same  as  in  the  previous  section 
Consider  each  different  condition. 


1.  Gross  weight 

=  K,  x  (48 

248,0001bs  240,0001bs 

=  2.0xl08  x(1.2)2’7 

=  3.27xl08 


x  (48,000/40,000)^ 


K  j  is  independent  of  changing  gross  weight. 

2.  Drag  index 

Kj  =  0.0135  +  0.000132  x  Drag  Index 
=  0.0135  +  0.000132  x  27 

-  0.01706 

K2  is  independent  of  changing  Drag  Index. 

3.  Altitude 

K1  =  K1DI27  x(pn,000ft  '  Pssl)1  2 

=  0.01706  x  (0.00 17007/0.0023769) 1,2 

=  0.01142 

K2'  Ss, 0001b  /(pl  1,000ft 'M1'2 

-  3.27xl08/(0.00 17708  0.0023769) *'2 

-  4.8893xl08 

Thrust  required  *  0.01142V2  +■  4.8893xl08  V2 

-  0.01142  x  (455xl.69)2  +  4.8893xl08  (455x1  69)2 

-  7,572  lb 


Compute  the  factor  C: 


Military'  thrust  =  21,800  -  0.4  x  Altitude 
=  21,800  -  0.4x11.000 
=  17.400  1b 

Fuel  flow; 

F,  =  20,840  -  0.37  x  (Altitude  -  8.000) 

=  20.840  -  0.37  x  (11.000  -  8.000) 

=  19.730 

p2  =  (Fj  +  333.3)  -  22.1  x  Drag  Index 
=  20,063.3  -  22.1  x  27 
=  19,466.6 

F3  =  F->  -  0.002  x  (Gross  weight  -  40,000! 
«  19,466.6  -  0.002  x  (48.000  -  40.000) 

=  19.450,6 

C  =  Military  thrust  Fuel  flow 
*  P.400  19450.6 
-  0  8945^ 

Thus 

Fuel  flow  »  Thrust  required  C 

-  ?,5'2  0.89457 

-  8.464  lbs  hr 


Now  the  fuel  flow  for  given  conditions  will  be  computed  by  interpolating  the 
hand  book  data.  Figure  3.7  shows  the  whole  diagram  to  compute  the  desired  fuel  flow 
of  given  conditions 

Computations  for  the  steps  are 
b  F  I  -  (8873  -*’988)x(2"'  -  20)  (40  -  20)  +  '988  -  829"5 
2 >  F2  -  (10324  -  9223)x(2'  -  20)  (40  -  20)  +  9223  »  %08.35 

3»  F3  -  (455  -  440 )x( 9608  35  -  829'  "5)  (480  -  440)  *  829'  '5  ~  8'S9  2' 

4 1  G1  -  (7729  -  6962 )X( 2'  -  20)  (40  -  20)  *  6962  *  '230  45 
b  G2  -  (8840-  '933)x(2' -20)  (40  -  20)  +  '933  -  8250  45 
6 1  G 3  -  (455  -  440 )x( 8250  45  -  '230  45X480  -  440)  *  '230  '5  -  '6 1  3  25 

■)  I  4  -  <  1 1000  -  <  non )xpH  3  25  -  8'89  23)  (  12000  -  8000)  *  8'89  2"'  =  '90'  24 

8.  HI  -  .92^9  .  8366ix(2'  -  20X40  -  2<>)  *  8366  -  867155 


DIAGRAM  FOR  DESIRED  FUEL  FLOW 


GROSS  WEIGHT  40,000  LBS 


'UDE  8,000  FT 


DRAG  INDEX 


0  |  7988  Fi  8873 


ALTITUDE  i: 

! ,  000 

FT 

A/S 

DRAG 

;  INDEX 

KTAS 

20 

27 

40 

440 

6962 

Gl 

7729 

80  9223  F2  110324 


480  7933  G->  8840 


|  GROSS 

WEIGHT  50,000  LBS 

ALTITUDE  8,000  FT 

11,000  FT  1 

ALTITUDE  12,000  FT 

KTAS  20 


440  8366 


DRAG  INDEX 


27  40 


Hx  9239 


55 


480  9522 


H2  10688 


A/S  DRAG  INDEX 

KTAS  20  I  27  4 0 


440  7447  Jx  8192 


455 


480  8315  J2  9216 


G/WT 

40,000  LBS 

A/S 

DRAG  INDEX 

KTAS 

27 

G/WT 

50,000  LBS 

A/S 

DRAG  INDEX 

KTAS 

27 

II 

55 

K4 

Figure  3.7  Diagram  for  Desired  Fuel  Flow. 
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Desired  fuel  flow  tor  given  condition  .s  s.242  'hs  hr 
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S.4M  -  S.242  i  S.242  x  1 1 W)  =  2  "  "  . 


I  v  (  DM  I  I  S|ON 


\  <  M*S  i  ,  ‘  \  *  1  ,1  \  \  *  % 


■*  r 


'  *’  .1*  .a.  f  v  »h*  '-vv  !• 


A  '  »  1  ' 


!'»■  *•*•>•* 
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•Jute'  'he  pr>.Kfhut\  i  making  i  mistake  r,  Merp«-,.i*i<in\  step* 

I;  ,s  recommended  fh.it  .i  similar  computer  program  Kc  dcseh.-pe 
ghi  stages,  i  e  .  fjWc-urt.  climb.  or  landing.  jsu.g  'he  same  concepts  ’ha 
:s  program 


“tV>v 


\PPKM>I\  A 

(  OMPI  1 1  R  PR(M,R  \M  (B  \SK  ) 


00310  */»/  9E*08  •  RfM  •  •  BASED  ON  5  W  AND  GROSS  WEIGHT  40000LBS  ** 
-ftrj’O  PRINT  -INPUT  A  i  T  I  T 1 1  Of  uROSS  WEIGHT. DRAG  INDEX, AIR  SPEED. A/S  INC* 

'ten  •  0  1  J  T  A  H  \  *  X  1  vP 

#i  -o  oils*  ooon.,»m  •  hem  >■  * i  varies  with  drag  index  ** 

1iO‘Q  0.  1  A  8'5L  6«A  ■  REM  ■  •  RAO  SCALE  BASE  •» 

9  r  0 » 0  rl  Oil  68S»D4  ’  UfM  »«  ABSOLUTE  TEMP  (R)  BASED  ON  ALTITUDE  ** 

,'.0  0  D5  C4**5  /Sbl  '  Rf,4  ••  PRESSURE  RATIO  BASED  ON  AITITUDE  »• 

ror-.RO  AO:(l  <.  •  |  -  ]  4  «  ’  •  T  1  )  •  REM  ••  SPEED  OF  SOUND  (SQUARED)  «• 

I0'<0  rviCN  D4  «»1  /  •  REM  ■  •  AITITUDE  CORRECTION  FOR  1 Q  •  r« 

M)  i  0  )  «i  «  l  •  *  9  '  REM  •«  P  A  ■  M  •  •  /  AITITUDE  CORRECTION  ** 

-Olio  i,'  «/*9  •  «EM  ■■  PA«p<«»/  ALTITUDE  CORRECTION  *• 

j  0  I  '  0  »,  «.•(«  <.0000  >**2  1  ' R  C  M*  *  WEIGHT  OTHER  THAN  40000LB  CORRECTION** 

0  0  1  *,  0  »  R  I  N  T  •  x  1  •  •  ,  *  1  •  x  /  -  '  .  *  2 

O-M’O  PRINT  TAB(  1  i  .  'Vfl  <  X  IS)'  .  TAB<  l  2  >  ;  'MACH';  T  A  B ( 2  0  )  .  'PAR  DRAG'; 

HJteO  PRINT  T  AB(  i<*  >  •  'PRO  DRAG' .  TAB(<*8  )  .  *MACM  DR  AG '  .  T  AB  (  6  2  >  ;  •  FUEL  FLOW' 

o  0  1  o  o  v  «’  v  1  ■  1  6  8  8  9  4 
0  0.’  0  0  M  v  2  '  A  9  •  •  0  5 

00.-1  3  08  1  '  REM  ••  C  DO  MACH  CORRECTION  BASE  »» 

03.20  IF  M >  8  THEN  48  0 

00/ JO  D 1  -  *  1 *v2**2  *  REM  •*  PARASITE  DRAG  (SUBSONIC)  ** 

00/40  D  J ’ D1 • D8  ’  REM  **  PARASITE  DRAG  (MACH  CDO)  ** 

00/50  D2 ’ K  2 / V2**2  ’  REM  »*  INDUCED  DRAG  ** 

OO/bO  DO  •  D  5 ♦ 02  *  REM  MACH  DRAG  ** 

OO/'O  D  7  -  D 1 ♦ D2  •  REM  »•  SUBSON IC  TOTAL  DRAG  »" 

00/80  T  21800  A*0  4  -  REM  •»  MIL  THRUST  OTHER  THAN  S.S.l.  CORRECTION  ** 
00/00  IF  A' -8000  THEN  S20 

00100  F  1  20840  -( A  8000 >*0  5  7  *  REM  F/FLOW  CORRECTT I  ON  ABOVE  8000  FT** 
00110  GO  TO  530 

305/0  FI *260000  64000*A  ’  REM  *«  E/FLOW  CORRECTION  BELOW  8000  FT  ** 

00  5  10  IF  1 9  ->20  THEN  360 

00  140  F  2 3  F  1  -  1  6  65*19  ’  REM  **  F/FLOW  CORRECTION  BELOW  DRAG  INDEX  20  ** 

00550  GO  TO  J7 D 

01540  F/=(F1*?55  51-22  10*19  !  *  F/FLOW  CORRECTION  ABOVE  DRAG  INDEX  20  * 
00570  IF  W<=41000THEN  400 

00580  F 5-F2-(W-40000)*0  002  !  REM  **  F/FLOW  CORRECTION  WITH  WEIGHT  »* 
00390  GO  TO  410 
00400  F5=F2 

00410  C 1 T/ F  5  1  REM  ««  COEFFICIENT  FOR  FUEL  FLOW  ** 

004/0  F  F  =  DO/C  !  REM  **  FUEL  FLOW  ** 

00450  PRINT  TAB(3);V1;TAB(8);M;TAB(18);D2;TABC32);D3;TAB(46);D0; 

00440  PRINT  TAB(60);FF 
00450  V1=V1*VP 
00460  IF  V l >6 1 0  THEN  510 
00470  GOTO  190 

00480  D9* . 55*<SIN((M-.8)/.8)»*3)/.0180 

00490  D8*1>D9 

00500  GOTO  230 

00510  PRINT  'C* ’ ;C. 'T*';T 

00520  STOP 

00530  END 


APPENDIX  B 

COMPUTER  PROGRAM  (HAND-HELD  CALCULATOR) 


*FF • 

51  P  L 

u:  pc.  ;*■ 

! 

52  • 

u:  f,l  ;■ 

0*  ’>1  <4 

5:  ST?  04 

u7  ♦ 

>*4  P[,  00 

54  PCL  04 

104  CTO  ir 

(■' 

5r  F'L  0T 

105  FCl  00 

0r  PPflHPT 

V  - 

10*  .4 

07  S'O  08 

57  STO  04 

107  • 

0 3  FCL  01 

53  2.8  E03 

10-:  chs 

0>. 

54  PCl  04 

104  21800 

10  PPO0FT 

t>0  f 

110  ♦ 

11  370  61 

61  STO  18 

Ill  STO  18 

12  PCL  *2 

62  PCL  01 

112  3000 

is  •»»»>• 

63  46000 

113  PCL  60 

14  PSOHPT 

64  / 

114  X.=  V 

15  STO  62 

65  2.7 

115  CTO  *PP* 

U*L6i  *W* 

6F  YtX 

lie  FCL  eo-- 

17  PCL  63 

67  PCL  10 

117  80*0 

tS  ‘VnATS:*’* 

6S  » 

118  - 

19  PFOHF'T 

69  STO  10 

119  .3? 

28  STO  83 

70  PCL  03 

120  * 

21  PCL  e2 

71  1.68394 

12!  CHS 

22  .8*0132 

72  * 

122  20840 

23  ♦ 

73  STO  11 

123  ♦ 

24  .8135 

74  PCL  87 

124  STO  19 

25  ♦ 

75  S6PT 

125  GTO  *11* 

26  STO  84 

76  t/X 

12(416L  ‘PP* 

27  6.375  f-06 

77  PCL  11 

127  PCL  80 

22  PCI  80 

78  * 

12S  .64 

2?  « 

7?  STO  12 

129  • 

36  CHS 

80  .8 

13*  CHS 

31  1 

31  V-V 

13!  260*0 

32  + 

82  GTO  *««• 

132  4 

33  STO  85 

834L8L  *RP* 

133  STO  19 

34  513.65 

84  PCL  11 

1344LEL  "LI* 

35  * 

85  Xt2 

135  28 

36  STO  06 

86  PCL  04 

136  PCl  02 

37  1716.5 

87  * 

137  X>Y? 

33  * 

88  STO  13 

138  GTO  *PD* 

39  1.4 

89  PCL  24 

139  16.65 

48  « 

90  * 

148  » 

41  STO  07 

91  STO  14 

Ml  CHS 

42  PCL  05 

92  PCL  11 

142  RCL  19 

43  EHTEPt 

93  XT2 

143  4 

44  5.2561 

94  1/X 

144  STO  20 

45  YtX 

95  RCL  10 

145  GTO  *EE* 

46  STO  *3 

56  * 

14C.4LBI  *DD* 

47  1116 

97  STO  15 

147  22.1 

43  Xt2 

98  PCL  14 

MS  * 

49  PCL  07 

99  ♦ 

149  CHS 

58  / 

100  STO  16 

158  333.3 

:c:  pa  u 

♦ 

1/4  ‘.TO  20 
!c-5n£i.  -ft' 
1  So  FCl  ¥  1 

157 

158  - 
!!4  .002 
Imp  • 

161  CHS 

162  FCL  :« 
US  ♦ 

164  STO  21 

165  1  X 

166  PCL  ie 
16?  * 

163  STO  22 
1(3  IX 
170  FCL  16 
1?1  • 

172  STO  21 

173  PCL  63 
17<  FIX  0 

175  -VUTSU 

176  APCl  X 

177  AVI  EH 
176  STOP 
173  PCL  12 
IS*  FIX  3 

181  *KACHr* 

182  APCL  X 

183  AVI  EH 

184  STOP 

185  PCL  15 

186  FIX.  1 

187  *C'I>=* 
183  APCL  X 
18?  AVI EH 
138  STOP 

191  PCL  14 

192  *D<P>=* 

193  APCL  X 

194  AVI EH 

195  STOP 

196  PCL  16 
19?  -D<n>=- 

198  APCL  X 

199  AVIEH 
208  STOP 


281  PC.  a 

2*>;  -fa* 

U*'  > 

2**»  u t  1  £  m 

?4C  •  T i’.P 
20*  *HtH  V*  ! 

2f*  PFOin 
28-:  x  --*1 
20H  :  rr.p 
2W  cto  *vv 
2!  1  *1 6L 

21 2  PCL  12 

213  .3 

214  - 

.8 

'16  / 

21?  EHTEPt 

215  3 

219  rtx 

220  S1H 

221  .55 

222  » 

223  .818 

224  / 

2C5  1 
2a  v 

227  SCO  24 
223  CTO  "FF* 
22?  .END. 
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